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ABSTRACT 



Context. The control of systematic effects when measuring background galaxy shapes is one of the main challenges for cosmic 
shear analyses. 

Aims. Study the fundamental limitations on shear accuracy due to the measurement of the point spread function (PSF) from 
the finite number of stars that are available. We translate the accuracy required for cosmological parameter estimation to the 
minimum number of stars over which the PSF must be calibrated. 

Methods. We characterise the error made in the shear arising from errors on the PSF. We consider different PSF models, from 
a simple elliptical gaussian to various shapelet par ametrisat ions. First we derive our results analytically in the case of infinitely 
small pixels (i.e. infinitely high resolution), then image simulations are used to validate these results and investigate the effect 
of finite pixel size in the case of the elliptical gaussian PSF. 

Results. Our results are expressed in terms of the minimum number of stars required to calibrate the PSF in order to ensure 
that systematic errors are smaller than statistical errors when estimating the cosmological parameters. On scales smaller than 
the area containing this minimum number of stars, there is not enough information to model the PSF. This means that these 
small scales should not be used to constrain cosmology unless the instrument and the observing strategy are optimised to make 
this variability extremely small. The minimum number of stars varies with the square of the star Signal-to-Noise Ratio, with 
the complexity of the PSF and with the pixel size. In the case of an elliptical gaussian PSF and in the absence of dithering, 2 
pixels per PSF full width at half maximum (FWHM) implies a 20% increase of the minimum number of stars compared to the 
ideal case of infinitely small pixels; 0.9 pixels per PSF FWHM implies a factor 100 increase. 

Conclusions. In the case of a good resolution and a typical Signal-to-Noise Ratio distribution of stars, we find that current 
surveys need the PSF to be calibrated over a few stars, which may explain residual systematics on scales smaller than a few 
arcmins. Future all-sky cosmic shear surveys require the PSF to be calibrated over a region containing about 50 stars. Due to 
the simplicity of our models these results should be interpreted as optimistic and therefore provide a measure of a systematic 
'floor' intrinsic to shape measurements. 

Key words. Gravitational lensing - Cosmology: dark matter - Cosmology: cosmological parameters 



1. Introduction 

Gravitational lensing by large scale structure (or 'cos- 
mic shear') has grown rapidly as a research field over the 
last d ecade (|Refregierll2QQ3bl : lHoekstrall2QQ3l : lMunshi et alJ 
20061 ). In the coming years, we can expect that the 
field will continue to mak e important contributions 



to cosmology. For instance, Albrecht et al. 
Peacock fc Schneid er (2006) have singled out cosmic shear 



( 2006h and 



as potentially the most powerful probe for constraining 
dark energy. 

A great deal of work is currently under way to develop 
techniques that will allow the maximum possible potential 
to be reached and not limited by systematic measurement 
errors. This can be done by designing and building in- 



struments with weak lensing as the central science driver. 
This approach leads to a top-down mission design pro- 
cess which begins with science requirements, which are 
translated into technical requirements and then into an 
instrument design. This method for tackling the problem 
of high accuracy weak lensing measurements is fundamen- 
tally different to the bottom-up approach that is currently 
used. In the latter approach an existing telescope is used 
and the science team is assigned the task of extracting the 
maximum possible information from their surveys (until 
they hit the systematic limit of the instrument). 

The field of gravitational lensing lends itself most nat- 
urally to top-down design, since most (but not all) of 
the systematic errors are not astronomical but are as- 
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sociated with the instrument and the atmosphere. These 
can therefore be controlled through instrument design and 
an optimised survey strategy. For this reason, many of 
the ambitious future imaging surveys that are currently 
under development have placed weak lensing as primary 
science driver, including the Dark UNiverse Exploreiu 
(DUNE), the SuperNovae Acceleration Prob^ (SNAP), 
the Panoramic Survey Telescope & Rapid Response 
Systenll (Pan-STARRS), the Dark Energy Surve^H (DES) 
and the Large Synoptic Survey Telescop^ (LSST). The 
issue is then to establish the instrumental requirements 
neede d to reach the full st a tistical potential of the survey. 

In lAmara fc Refregierl (|2QQ7af ). a wide range of sur- 
vey parameters, such as area and depth, were considered. 
Their impact on the statistical potential of a lensing sur- 
vey was calculated and summarised in a scaling relation 
that can be used for survey designs to trade-off one prop- 
erty of a survey against another. In brief, this work finds 
that, once the median redshift of a survey is sufficiently 
high {z > 0.7), the optimal survey strategy is to make the 
lensing survey as wide as po s sible. Similar results have 
been found bv iHeavens et~al In the same spirit, 

Amara fc Refregierl ( 2QQ7bh looked at the requirements an 
ultra- wide field survey places on the control of systematics 
and concluded with a set of scaling relations that show the 
tolerance on residual systematic errors as a function of sur- 
vey parameters. In particular, two types of shape measure- 
ment systematics were considered: multiplicative and ad- 
ditive systematics which are, respec tively, correlat e d and 
uncorrelated with the lensing signal. !Huterer et al ] (|2006) 
have also studied the impact of multiplicative and additive 
errors. 

In the present study, we link the above systematic re- 
quirements to errors associated with measurement of the 
point spread function (PSF) of the instrument. Indeed, 
since galaxy images need to be PSF-corrected before 
shapes can be measured, errors in the estimation of the 
PSF are propagated into errors on galaxy shapes and can 
mimic the shear signal. The PSF calibration is driven by 
3 factors: 

1. The PSF model: The model that is chosen to de- 
scribe the PSF will never perfectly describe the true 
PSF hence the choice of PSF model introduces an 'a 
priori' systematic; 

2. The interpolation scheme: The PSF will vary over 
the field which means that the PSF at a galaxy position 
needs to be interpolated from the PSFs of the nearby 
stars; 

3. The finite information available from each star: 

A star provides an image of the PSF that is noisy and 
pixelated. Thus, the PSF information that we are able 
to extract from each star is finite. As a result, to reach 



http:/ /www. dune- mission, net 

http://snap.lbl.gov/ 

http: / / pan-starrs.ifa.hawaii.edu 

http:/ /www. darkenergysurvey.org 

http:/ /www. lsst.org 



high precision, it is necessary to combine the informa- 
tion from several stars. 

This paper is devoted to the last factor. We quantify the 
accuracy of the PSF calibration with analytical predic- 
tions in the case of infinitely small pixels and use image 
simulations to quantify the pixelation effects. Previous 
work has also studied the impact that PSF errors are 
likely to indu ce in cos n iic sh ear measurements. For in- 
stance, Stabe nau et al.l (|2QQ7l ) looked specifically at the 
PSF from the current SNAP design and translated this 
into systematic errors. The study that we present here 
sets out to be more general and uses simplified PSFs to 
try and quantify the systematic floor of an instrument. 
We investigate the extent to which the pixel scale will de- 
grade the amount of information relative to our analytic 
predictions but we do not explore the optimal methods for 
co mbining multiple e xposures, su ch as the one proposed 
bv I Jain et a D (l2006l). iHigh et all ([2007) have studied the 
impact of the pixel scale on the number of useful galaxies 
(assuming the PSF is known perfectly). While they find 
a weak dependence, we expect that the driving factor for 
pixel scale is shape measurement systematics rather than 
statistical errors. 

For a sur vey with a given statistical potential we use 
the results of Amara fc Refregier (|2007bh to find the upper 
limit on the systematic errors in the shears such that the 
induced bias in cosmological parameters remain subdomi- 
nant compared to the statistical (marginalised Fisher ma- 
trix) errors when estimating cosmological parameters. We 
then convert the requirements on shear systematics into 
the minimum number of stars (A^*) needed to measure the 
PSF to the level of accuracy required. On scales smaller 
than the area containing TV*, there is not enough informa- 
tion from the stars to calibrate the PSF. Therefore, the 
systematic errors in the cosmological parameters can be 
dominant over the statistical ones. This means that these 
small scales should not be used to constrain cosmology un- 
less the PSF variability is known to be extremely small. 
clearly depends on the Signal-to-Noise Ratio distribution 
of stars, since bright ones contain more information about 
the PSF than faint noisy ones. It also depends on the sta- 
bility and the complexity of the PSF, characterised by the 
number of degrees of freedom that must be estimated or 
interpolated from the stars. Each degree of freedom in- 
creases N^. 

Figure [T] shows a graphical example to illustrate the 
calibration requirement we are considering. We see a cen- 
tral galaxy surrounded by several stars. The grey shaded 
region shows the area over which the PSF must be cali- 
brated and in this example contains 11 stars. Stars do not 
provide enough information to model the PSF variations 
on scales smaller than this area. 

This paper is organised as follows. Section [2] explains 
the weak lensing and cosmological context and summarises 
the issue of accurate PSF calibration. Section [3] shows the 
analytical predictions of the PSF calibration accuracy, ex- 
pected in the case of infinitely small pixels. Section [H de- 
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Fig. 1. Illustration of the required number of stars 
characterised in this paper. When measuring the shape of 
a galaxy (marked by a red spiral), the PSF needs to be 
calibrated with at least nearby stars (black asterisks) 
contained in the shaded region. In this example = 11. 
On scales smaller than this there is not enough informa- 
tion coming from the stars to measure the PSF variations. 



scribes the simulations we use to validate these predictions 
and extend them to finite pixels. In section [5] we give the 
final accuracy of the PSF calibration and derive the num- 
ber of stars it requires. 



2. Weak lensing 

2.1. Shear Measurement 

Weak gravitational shear is locally estimated using the 
shapes of background galaxie s (for a review, see e.g. 
Bartelmann fc Schneider 2QQlh . Shear estimation meth- 
ods can be divided into two families: those comput- 
ing an estimator of the shear from a set of weighted 
sums over pixel values, for instance the common 
J£SB+' method (Kaiser e t al. 19951 : iLuppino &: Kaiser 
19971 : iHoekstra et al. 1998) . and those fitting a model to 
the observed galaxy shape and deriving an estima tor from 
the m o del. This includes the meth o ds proposed by Kuiiken 
(Il999h: iBernstein fc Jar^ (|2QQ2h : iRefregier et al.. (.2002 ) 
and lVoigt fc Bridle! tOoK 

In this paper, we need a general formalism to prop- 
agate the error on the PSF into an error on the shear 
estimate. For this purpose, we consider as shape parame- 
ters the 2 component ellipticity e and the squared radius 
i?^, both defined using the unweighted second order mo- 
ments of the galaxy. For an object with surface brightness 
/(xi,X2), the total flux F^^^ is the zeroth order moment 
of the surface brightness: 



F(°) = I d'xfix), 



(1) 



the centroid x^^^ is given by the first order moments di- 
vided by the total flux: 



(1) 



^ J (fxxif{x) 



(2) 



and the quadrupole moment matrix is given by the second 
order moments divided by the total flux: 



(2) 



1 



d^X {Xi - X^^^''){Xj - xf"") f{x) . 



(3) 

The square rms radius and the 2 component ellipticity 
€ = [ei, 62] are defined by: 



ei = 



Qii 
Qu 



C22 
^22 



^22 



2Q 



12 



Qii + Q22 



(4) 
(5) 

(6) 



Consider now we have unbiased estimators of E? and e 
with variations 5E? and 8e around the true values. We 
adopt the foUowing definitions: 



(7) 
(8) 



The factor 2 in equation [8] is due to the fact that e has 2 
components and we define cr[epsF] as the standard devia- 
tion of one of the components. 

The weak gravitational shear 7 = [71, 72] can then be 
shown to be estimated using: 



7 



(p7)- 



where 



= 2-<|€gair >^1.84. 



(9) 



(10) 



The symbol ^' indicates an estimator so that < 7 >= 7, 
is the shear susceptibility and the subscript 'gal' cor- 
responds to values measured on a galaxy. The value of 



< l^gail > ^ 0-4 comes from the typical ellipticity dis- 
tribution in current data sets. 

In practice, the measurement of the galaxy elliptic- 
ity is uncertain because of the noise in the image. We 
write the induced error as ^e"^^^®, which has a null aver- 
age < ^e"°^^® >= 0. There may also be a systematic effect 
Se^^^ so that the estimated ellipticity is 

egal = €gal + 5e^^^ + 5€"°"^ (H) 

In practice Se^^^ may be due to s everal factors. For in- 



stance, the STEP collaboration (iHevmans et al. 20061 : 



'Massev et aP I2QQ7I: iRhodes et aP l2008[ ) investigate the 
contribution to 5e^^^ from imperfections in the shear mea- 
surement methods. In a different spirit, this paper focusses 
on the contribution to 6e^^^ from the limited information 
available on the PSF due to the photon noise in the star 
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images. We identify the variance ( |(5e®-^®| ) with the quan- 



tity cTgyg examined in Amara fc Refregier ( 2QQ7bh . which 
indicates the variance of the systematic errors 

-2 /i^esys|2\ 



,2 _ 



(p7)- 



(12) 



We translate an upper limit on a^^^ into constraints on 



This equation confirms the intuition that the PSF elliptic- 
ity of a cosmic shear survey should be small. We see that 
the last ellipticity term inside the brackets is ^|6psf|^^ 
and should be reduced when optimising the survey to have 
epsF as small as possible. 

For instance, to reach the requirements of (jg < 10 



the number stars that are needed to cahbrate the PSF. T^*^ a ^typical well sampled cosmic shear survey with 

= 0.16, Rgai > l-5i?psF and epsF ^ 0.05, would 



2.2. System atics from PSF calibration 

Since a galaxy image needs to be PSF corrected before 
its shape can be measured, errors in the estimation of the 
PSF propagate into an error Se^^^ in the measured elliptic- 
ity. Here, we consider errors in the PSF radius parameter 
RpsF and in the 2 component PSF ellipticity €psf- To 
first order and for the unweighted moments (definitions 
given by equations [3] to [6]) , we have (see Appendix [Al for 
details): 



(Cgal — epsF. 



5 {Rpsf) [ RpsF 



RL 



Rcr 



6epsF (13) 



where Rg^i and i^psF are the radius parameters of the 
galaxy and the PSF respectively. The first two terms show 
that systematics due to an error on the PSF size are pro- 
portional to the ellipticities of the galaxy and of the PSF. 
In the following, we show that the latter should be op- 
timised to be as small as possible. The last term shows 
that systematics due to an error on the PSF ellipticity are 
proportional to the squared ratio between the PSF and 
galaxy sizes. 

Combining equations [12] and [13] gives the propagation 
of PSF errors into dgyg. For this purpose we make the 
following simplifications: 

1. The galaxy is not correlated with the PSF (i.e. the 
crossed-terms < €gai.€psF > and < €gai.^6psF > are 
equal to 0). 

2. The error on the PSF ellipticity (^€psf) and the PSF 
ellipticity itself (epsp) are not correlated. This is war- 
ranted by the fact that, in the assumed small ellipticity 
regime, ^€psf does not have any preferred direction, 
implying < epsF-^epsF >= 0. 

3. We assume that the ellipticity and the inverse squared 
radius of the galaxy are also uncorrelated. More ex- 
actly we assume: 

2\ 



Rt 



l/Rf 



gall 



(14) 



This assumption is reasonable for our work on the PSF 
calibration presented in this paper. 

With these simplifications, we obtain: 
f Rpsf\^\ 



2(T"[epsF] 



((!< 



€PSF 



a[R^ 



PSFJ 



RpSF 



(15) 



require 



^[-^psf] 

rtpgp 

c^[epsF] 



< 10-^ 

< 10-^ 



(16) 
(17) 



In the following, we discuss how these upper limits trans- 
late into requirements on the PSF calibration. 

3. Analytical model 

In this section, we consider the general problem of fitting 
a 2D model of the PSF to a pixelated image of a star with 
additive uncorrelated gaussian noise using minimisa- 
tion. This allows us to derive a number of analytic results 
in the limit of infinitely small pixels (i.e. infinitely high 
resolution), which will serve as a useful comparison base 
for the numerical simulation studies presented in section 

m 



3.1. General 2D fit 

Consider the PSF surface brightness as a function of po- 
sition on the image x to be described by a model m{x;p) 
parameterised by parameters p. The observed surface 
brightness of a star is f{x) = m{x] p(PSF)) + n{x) where 
p(PSF) is the true values of the parameters, n{x) is 
the noise, which is assumed to be uncorrelated (from 
pixel to pixel) and gaussian with < n{x) >= and 
< n{x)'^ >= cr^ is assumed constant across the image. The 
usual x^-functional is given by 

X'(p) =^cT-2[/(xfc)-m(xfe;p)]% (18) 



the sum being over all pixels k in the image. The usual 
estimator p is constructed by requiring that dx^/dp = 
when evaluated at p. The covariance matrix of the esti- 
mated parameters is given by the inverse of the Fisher 
matrix: 

coy[p,,Pj]^{:F-% (19) 



with 



2 dm{xk; p) dm{xk;p) 



dpi 



dpj 



(20) 



and the variance of a function P(p) of the fitted parame- 
ters is: 

-^[^]-E^a^-b..,]. (21) 
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In the following, we consider 2 such functions of the pa- 
rameters p: the rms radius squared and the 2 compo- 
nent ellipticity e (as defined in equations [H to (H) , in the 
case of a simple elliptical gaussian PSF (section [321) 
in the case of more complex PSFs defined with shapelets 
(section [XS]) . We define the associated dimensionless com- 
plexity factors iIjji2 and i/j^ such as: 



a[R' 



1 
S 



where S is the Signal-to-Noise Ratio defined as: 

i?(o) 



(22) 
(23) 

- - .[F(0)] • (24) 

cf[F^^^] is the standard deviation of the total flux. 
Basically, ?/;^2 and ipe characterise the numbers of degrees 
of freedom in the PSF model associated with R^ and e. 

In the limit of infinitely small pixels, the sum over pix- 
els in equations [TSl and [2Ql can be replaced by a continuous 
integral over the object. A number of analytical results can 
be derived. In section [321 we derive these for an elliptical 
gaussian PSF and in section [331 we study more complex 
PSFs described with shapelets. 

3.2. Elliptical gaussian 

We first consider a 2D elliptical gaussian model parame- 
terised as: 

A 



m{x; p) - 



exp 



1 



X 



^Q-^^x-x"") 



p = {x'^,A,ai,a2,(y) 



(^5) 
(26) 



where a\ and 02 are the rms major and minor axes of the 
gaussian, respectively, A is a parameter which controls the 
amplitude, is the (true) centroid, and stands for the 
transpose operator. The total flux (as defined in equation 
El) is: 

= Ay/^^ , (27) 
the centroid (as defined in equation |2]) is: 

rc^" = x1 (28) 
and the quadrupole moment (as defined in equation |31) is: 



Q = R(a)^ 



ai 
al 



R(a) 



(29) 



where a is the position angle of the major axis counter- 
clockwise from the x-axis and: 



R(a) = 



cos a sm a 
— sin a cos a 



(30) 



is the rotation matrix which aligns the coordinate system 
with the major axis. The ellipticity and the squared radius 
(as defined in equations |4| to |6l) are: 



i?2 



(31) 
(32) 



This parametrisation is particularly convenient be- 
cause for infinitely small pixels the Fisher matrix of the 
parameters, given by equation [20l is diagonal with diago- 
nal elements: 



1 



aY A^^ 0^2^ 



1 1 {a\-alf 



2a\a2 



(33) 



Consequently, with equation [ 
and minor axes are: 



the errors on the major 



S' 



From equations [31} [32l and [34l it follows that: 



cj[R^ 



^2 ' 



4v^ 



i?4 S 



(34) 

(35) 
(36) 



We generalise this simple elliptical gaussian model to more 
complex PSFs in the following section and test these equa- 
tions for finite pixels using simulations in section [H 

Under the approximation of small ellipticity (e < 0.1), 
the dimensionless complexity factors ?/^^2 and ijje (defined 
by equations [22l and [23l) are constant (i.e. do not depend 
on the object). Indeed, e < 0.1 implies ai^_a2_and there- 
fore in equations [35l and [36l we have Y^of + of :^ R^ /\/2 
and Aa\a2 — R^- This leads to: 



'0i?2 (gauss) = ?/^e (gauss) = V2 . 



(37) 



3.3. Shapelet model 



To explore a wide variety of possible PSFs we consider a 
basis set that allows for complexity. In this framework, the 
shapelet basis sets are particularly convenient because: (i) 
shapelets provide some orthonormal basis sets that have 
already been studied in num b er of p u blications (.Refregier 



tpi 

2003al : [Refregier fc BaconI [20031 : [Massev fc R.efregier 
2005[ ) tested on the STEP2 and STEP3 simulated 
data ( Massev et al. I l2nn7l : [Rhodes et a].ll2nn8h and used 
for the weak lensing analysis of the Canada-France- 
Hawai-Telescope-Legacy- Survey (CFHTLS) data set 
in the framework of a comparison with X-ray surveys 
(jBerge et al.[[2007| ): (ii) as we show in the following, in 
the case of simple objects the standard deviations of the 
squared radius cr[i?^] and ellipticity a[e] can be written 
as proportional to some complexity factors containing all 
the information about the basis. 

A shapelet basis set is characterised by 2 parameters: 
^max, the maximum order of the functions in the basis and 
the scale parameter p. The surface brightness f{x) of an 
object is described by: 



(38) 



n=0 m= 



with Xn,m the polar shapelet functions and fn,m the (com- 
plex) polar shapelet coefficients of the object. The scale of 
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the oscillations described by the function Xn,m is propor- 
tional to l/(n + |m|). It is often convenient to impose a 
lower limit to the scales which are described. This corre- 
sponds to setting the coefficients with |m| > nmax — ^ to 0. 
This configuration is called 'diamond' and is only defined 
for an even rimax- For simple objects with reasonable sub- 
structures and tails, the number of coefficients required is 
small when ~ 2/?^. This relation is exact for a circular 
gaussian represented with n^iax = 0. 

We show in Appendix [B] that, under the approximation 
of small ellipticity (e < 0.1) and for stars described with 
a judicious scale parameter where = 2/3^, the dimen- 
sionless complexity factors defined by equations [22l and [23l 
{cr[R'^] = R'^iIjji2/S and a[e] =iIj^/S) depend only on the 
basis (i.e. not on the object itself) and are given by: 



^R^ ^ ]Jj{N^l), (39) 



V^e - y— (A^ + 4) without diamond , (40) 



V^e - y—{N-2) with diamond (41) 

where is the largest even integer lower than or equal to 
^max- Figure [3] illustrates equations [22l [23] and [39] to [41] 
by showing (j[e] and cr[R'^] as a function of nmax- V^e and 
ipji2 are given in table [2] for 16 bases with nmax between 4 
and 16 and for the elliptical gaussian model presented in 
the previous section. 



4.1. Elliptical gaussian PSF 

Figure [2] shows the measured standard deviations on R^ 
and e as a function of <S, for a range of pixel scales (from 
0.9 to 3.9 pixels per PSF FWHM). The black diamond 
symbols correpond to the highest resolution simulations 
(with 3.9 pixels per PSF FWHM). These points are close 
to the predictions (black curve) which are calculated an- 
alytically for the idealised case with infinitly small pixels 
(i.e. inifinitly high resolution, see section [3^ . The upper 
and bottom panels show different representations of the 
same curves: on upper panels the y-axis shows a, while on 
the bottom panels, it shows the relative differences with 
the infinite resolution case. As the pixel size increases (i.e. 
when the resolution decreases), the standard deviation in- 
creases monotonically and moves away from the analytical 
predictions. This shows that the analytic model can de- 
scribe a high resolution case while simulations are required 
to take into account pixelation effects. 

Our results in this paper are based on analytical pre- 
dictions and assume the optimistic case where pixelation 
is good enough not to degrade the results. From Figure [2] 
we see that, for the smallest pixel scale we consider (3.9 
pixels per PSF FWHM), the results from the analytic cal- 
culations and the simulations agree to better than 10%. 
As the size of the pixels is increased to 0.9 pixels per PSF 
FWHM we see a dramatic increase of the standard devi- 
ation of a factor ~ 10. Some of this pixel scale effect will 
be mitigated through the use of dithering, but more sim- 
ulations are needed to fully quantify the extent to which 
this will help. 



4. Simulations 

In this section, we use image simulations to test and vali- 
date the analytic predictions of the previous section. These 
simulations allow us to investigate the effects of pixela- 
tion in the case of an elliptical gaussian PSF. In this sim- 
ple case we show that pixelation degrades the accuracy of 
the shape measurements and thus, if the pixel scale is in- 
creased, the number of stars needed to calibrate the PSF 
also increases. 

Star images are simulated for: (i) a simple elliptical 
gaussian PSF and (ii) complex PSFs der ived from the 
STEP HI simulations ([Rhodes et al.ll2008[ ). In ah cases 
the pixelated images are produced by distributing the star 
centroids randomly and uniformly over the central pixel. 
Additive gaussian noise is then added to each image and 
the model fitting is done using minimisation (for the 
work presented here, we have not added Poisson noise, 
which arises from the object itself), gaussian stars are sim- 
ulated at high resolution (each pixel is considered as the 
sum of 7 X 7 sub-pixels), while shapelet stars are simulated 
using the analytic pixel i ntegration scheme proposed in the 
online IDL pipelin^ bv lMassev fc Refregier ( 2005 ). 



http:/ /www. astro. caltech.edu/^ rjm/shapelets 



4.2. Shapelet PSF 

We now turn our attention to the expected standard de- 
viations of R^ and e (predicted by equations [22] and [23|) 
for more complicated PSFs. We still work under the setup 
that the model fitted to the data can exactly describe the 
truth. With the notation used in section 13.11 this means 
that when we fit the parameters p of a model m(ic,p) 
there exists a solution p(PSF) that exactly describes the 
PSF. 

We model the PSF by decomposing PSF-D of STEP 
III (|R.hodes et al.ll2008h into shapelets. We create 16 PSF 
models with increasing levels of complexity by decompos- 
ing the STEP PSF into basis sets with increasing nmax 
values. The fitted models are listed in table [2l This gives 
16 slightly different PSFs close to PSF-D (e^ vary between 
-0.02 and 0.02) which are shown in figure [4] 

The effect of the level of complexity of the PSF model 
on the standard deviation of a measurement of the size 
and ellipticity of the PSF is found by generating ~ 10^ 
realisations in the interval 10 < <S < 1000 for each of 
the 16 PSF shapelet basis sets and by fitting with the 
corresponding shapelet model. 

We compare the analytical predictions given by equa- 
tions [22l [23] and [39] to [4l] with simulations in figure [3] We 
see they provide a good description, although they some- 
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X 0.9 pix. per PSF FWHM 

□ 1.1 pix. per PSF FWHM 

A 1 .4 pix. per PSF FWHM 

+ 1 F pi-. . p'-r P-F FWHM 

+ 2.3 pix. per PSF FWHM 

O 3.9 pix. per PSF FWHM 




SNR SNR 



Fig. 2. Solid lines: analytical predictions in the case of infinitly small pixels for a[e] and cr[R'^] (see equations [35l 
and [36]). Data points joined by straight lines: simulation with finite pixel scales. All these results are for an elliptical 
gaussian PSF with an ellipticity e :^ 0.1. Each data point corresponds to about 10^ realisations. Black diamonds, red 
horizontal-vertical crosses, orange stars, green triangles, blue squares and purple diagonal crosses correspond to sizes 
ai = 1.7, 1, 0.7, 0.6, 0.5 and 0.4 pixels respectively. This corresponds to about 3.9, 2.3, 1.6, 1.4, 1.2 and 0.9 pixels 
respectively per PSF FWHM. Each data point corresponds to about 10^ realisations. On the upper panels the y-axis 
shows a[R^]/R^ and (Je , while on the bottom panels it shows the relative differences as compared to the infinite 
resolution case. The Signal-to-Noise Ratio is noted 'SNR' rather than S in the text. 



times underestimate (j[R'^] and a[e] by a few percent. The 
results are similar whether or not the diamond configura- 
tion is used. 



5. Requirements 

In this section, we estimate the number of stars re- 
quired to estimate the PSF to sufficient accuracy for a 
given survey. First, each star is an independent realisa- 
tion of the PSF. Thus the PSF size and ellipticity can be 
estimated by: 



i?2 



PSF 



epsF 



Ek^^'lRl] EkS! 
T,k ^k/(T'^[ek] ^ Efc gfc si 



(42) 
(43) 



where the sums are made on all the stars. These estimators 
are the minimum variance weighted averages if we assume 
the probability distribution in each of R'^ and the two 
components of e is a gaussian for each individual star. 
Second, variances are given by equations [22] and [23] and 
can be written: 



i?2 



(44) 
(45) 



We have defined an effective signal-to-noise of the stars 
5eff as the expected rms signal-to-noise given the max- 
imum and minimum limiting values of S used for PSF 
calibration (<Smax and «Smin) and the number density of 
stars per unit of S: dn/dS^ in the data set: 



s? 



'eff 




dn 
dS 



S'dS 



We then define the effective number of stars A^* by 



(46) 



(47) 



where the sum is over the stars available around the 
galaxy. Therefore, the effective number of stars is equal 
to the actual number of stars if all the stars have the same 
S. 

Substituting into equation [15] gives the requirement 
on the number of stars needed to calibrate the PSF. If 
the PSF ellipticity is at most a few percent then we can 
neglect the ( |6psf|^ / term leaving: 



7 cr^^m\ 



sys> 



^eff 



fRpsFV 



(48) 
. The num- 
ber of stars is a central requirement of cosmic shear 



where a^^^ is the upper limit acceptable for a^ys 
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Fig. 3. Lines: analytical predictions in the case of infinitely smah pixels for (j[R^] and cr[e] (see equations [22l and 
[23|) . Data points: simulation with small pixels (i.e. 10 pixels per PSF FWHM). All these results are for stars imaged 
through the complex PSFs shown in figure [4] and with 10 < 5 < 1000. Each data point corresponds to about 1.5 x 10^ 
realisations. Triangles and diamonds correspond to basis sets with and without the diamond configuration respectively. 
For cf[R^] the solid line corresponds to the theoretical prediction (which does depend on the diamond option). For 
cr[e], solid and dashed lines correspond to theoretical predictions with and without the diamond option respectively. 
Without the diamond option, the complexity factors depend on the largest even integer lower than or equal to nmax 
(see equations [39l and [4Q|) . this is why theoretical predictions look like a 'staircase'. On the other hand this behaviour 
does not appear with the diamond configuration because the latter is defined only for even values of nmax- We see 
that analytical predictions provide a good description, although they sometimes underestimate cf[R^] and (j[e] by a 
few percent. The Signal-to- Noise Ratio is noted 'SNR' rather than S in the text. 



surveys: if we need to measure the PSF so that the sys- 
tematic effects of the shear power spectrum stay below 
cTg^^-, we need to combine the information from at least 

stars. Thus, for a given survey, one needs to estimate 

in order to optimise the instrument and mission de- 
sign. This is a requirement on the PSF stability. Indeed 
the star density must be taken into account to give the 
minimum area containing N^^ stars, over which the PSF 
must be stable. 

Low S stars do not contribute a great deal of in- 
formation to the PSF calibration. For instance if the 
density of stars scales roughly as dn/dS{S) ex 1/5, 
w hich is con sistent with the simple star count model 
of iBahcall fc Soneiral (|l980l ) , equation [46] shows that the 
bright stars strongly dominate. In fact, not only do low 
S stars bring little information, they may also induce a 
bias in the PSF calibration. This bias will be studied in 
forthcoming work. For the moment, we avoid this regime 
by adopting an arbitrary lower limit of 5min = 100. 

On the other hand, the high S cut off 5max and the 
star density distribution dn/dS^ depend on the proper- 
ties of the surveys. They depend on a number of factors 



including: the line of sight in the Milky Way, the instru- 
mental configuration and the observing strategy. Equation 
[48l therefore, can not be computed for a general case, but 
it can be simplified and scaled to typical values, to give 
easily readable requirements on the star population in a 
data set. First, since the galaxy size distribution is steep, 
it is pessimistic but reasonable to approximate: 

((ls)>[(ls)J' 

where (i^gai/^PSF)min is the minimum value that this ra- 
tio can reach, typically about 1.5. Second, we note that the 

variance ^|€gai|^^ is around 0.16, as stated in section [2j 
thus for usual PSF models where the complexity factor of 
the squared size ?/^^2 is of the same order as the complexity 
factor of the ellipticity ?/^e the expression inside the brack- 
ets in equation [48] is driven by ^jj^ and we can neglect V^^2 . 
For example, for non-diamond shapelets with nmax = 4 or 
for diamond shapelets with nmax = 6, the ellipticity com- 
plexity factor is 3 (see table [2]). Third, as mentioned 
in section [2| .Amara fc Refregier (2007bi) find for a DUNE- 
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like cosmic shear experiment (cTg^^ )^ = 10~^. Fourth, we 
take for 5eff a typical value of 500, which is roughly S of 
an AB magnitude 20 star when imaged with a DUNE-like 
telescope (1500 seconds in a broad RIZ band). With these 
central values, we obtain: 



(-figal/-^PSF)n 

1.5 



10- 



(50) 

Therefore the PSF of a DUNE-like survey will need to be 
stable over a region containing about 50 stars. 

The scaling relation given by equation [50] allows one 
to study different survey configurations. This is illustrated 
by table [1] that gives A^* for 6 typical configurations com- 
parable to: the Canada- France-Hawai- Telescope Legacy 
Survey (CFHTLS), the Kilo-Degree Surve}0 (KIDS), the 
Large Synoptic Survey Telescope (LSST), the SuperNovae 
Acceleration Probe (SNAP) and the Dark UNiverse 
Explorer (DUNE). To build this table, we have estimated 
the accuracy a[wo] and a[wa] that can be achieved ac- 
cording to the surface covered, the median redshift Zm 
and the galactic surface density Ug in the case where sys- 
tematics are lower than statistical errors and when fit- 
ting the 7 cosmological parameters: Qrn^ crs, wq^ Wa 
and n. This estimation is made through the s caling rela- 
tions proposed by Amara &: Refregierl (|2007a[ ). We have 



then estimated the value of (cTg^^')^ required to achieve 
th is assumption of su b-dominant systematics, according 
to lAmara fc Refregierl (| 2 00 7b ) . Finally we have computed 
A^* according to our scaling relation (eq. [50]) . 

1. The two first lines show surveys comparable to the cur- 
rent largest data sets optimised for cosmic shear. For 
instance the C FHTLS which currently covers 50deg^ 
( Fu et al.ll2007[ ) and will eventually cover 170 deg^ The 
relatively poor constraints got on wq with such surveys 
(and the absence of constraint on Wa) illustrate the fact 
that cosmic shear surveys of the current generation 
do not constraint significantly the dark energy. They 
rather aim to constraint the combination as x 
However t^iis does not change tjie requirement on a^^^- 
since lAmara fc Refregierl (j2007b ) show that this re- 
quirement does not depend on the considered cosmo- 
logical parameter. One can see that these surveys need 
to calibrate their PSF over few stars. Assuming a typ- 
ical star density of ~ l.arcmin"^ at these magnitudes, 
this shows that at scales smaller than few arcmin the 
PSF correction may introduce signific ant systematics. 
This is consistent with the results of [Benjamin et al.l 
( 20071 ) who describe a joint analysis of 100 deg^ from 
several surveys and find significant B modes on these 
scales. 

2. The third line corresponds to a cosmic shear survey of 
the next generation, such as KIDS/VIKING. 

3. The fourth line shows an ambitious deep survey from 
the ground covering 15,000deg^, which is the largest 



cosmic shear survey that can be perfomed from Mauna 
Kea with a small airmasfl This illustrates one of the 
main limitation of doing cosmic shear from the ground: 
although the survey is rather deep, the galaxy sur- 
face density Ug and the median redshift Zm are lim- 
ited. This is due to the fact that faint galaxies are 
too small: they have sizes comparable to the PSF and 
can not be included in a cosmic shear analysis. This 
is included in the scaling relation [50] through the min- 
imum dilution factor (i^gai/^PSF)min- Note that it is 
not possible to change significantly the reference value 
of (i^gai/^PSF)min = 1-5 adopted in equation [50] be- 
cause it is at the power 4. For instance, a change from 
1.5 to 1.2 to increase Ug would also increase by 
a factor 2.5. Such a survey i s sini ilar to the LSST 
(|Kahn fc LSST Collaboration! liooi ) except that the 
LSST collaboration is developing a specific observing 
strategy that aims to drastically reduce i/j^. This strat- 
egy is based on stacking many very short exposures 
taken with different orientations of the camera in or- 
der to circularise the PSF. Another effect of this strat- 
egy is that the effective Signal-to-Noise Ration 5eff of 
stars may be much larger than 300, as assumed in ta- 
ble [TJ But this strategy also has a fundamental limit: 
the more t/je is reduced, the more ipji2 increases (for 
instance because of the astrometric errors while stack- 
ing the low Signal-to-Noise Ratio exposures). For such 
a LSST-like strategy, estimating A^* through equation 
[50] is not relevant because in this equation we assume 
that tjje and ?/^^2 are of the same order. must be 
estimated through equation [48] 

4. The fifth line corresponds to a deep space survey cov- 
ering 1000 deg^. This is simi lar to SNAP as described 
by Massev et all (|2004l ) and iR^efreder et all (|2004). 

5. The last line corresponds to a DUNE-like survey (full 
sky from space with a medium depth and a reasonable 
galaxy density) adopted as the reference in the scaling 
relation (eq. [5Q|) . 

6. Conclusions 

We have studied the PSF calibration requirements for cos- 
mic shear to measure cosmological parameters. We con- 
nect the finite information that we are able to extract from 
stars to the statistical error of the PSF calibration, then to 
the error on the galaxy ellipticity estimation. We express 
our results in terms of the minimum number of stars (A^*) 
that are needed to calibrate the PSF in order to keep the 
systematic errors below statistical errors for cosmological 



|http: / / www.eso.org/ sci / observing/policies / cSurveys / sciencecSurveys.html 



^ On one hand, half of the sky is available for cosmic shear, 
corresponding to ^20,000 deg^. The other half is masked by the 
Milky Way. This is why one uses to call 'full-sky' the surveys 
covering 20,000 deg^. On the other hand, in order to maintain 
the air mass in a reasonable range, we make the optimistic 
assumption that the minimum elevation of ground observations 
is 50 deg in local horizon coordinates. With this constraint, a 
telescope located at Mauna Kea (i.e. Hawaii) can cover about 
15, 000 deg^ during the year 
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surface 
(deg^) 




rig 

(.arcmin"^) 


a[wo] {cr[wa]) 






Ground surveys 

(5eff = 300) 


50 


0.8 


10 


0.9 (nr) 


4 X lO""" 


3 


170 


0.8 


10 


0.5 (nr) 


2 X 10"^ 


6 


1,500 


0.8 


20 


0.1 (0.6) 


6 X 10"^ 


25 


15,000 


0.9 


30 


0.03 (0.1) 


10-^ 


104 


Space surveys 
(<Seff = 500) 


1,000 


1.2 


50 


0.07 (0.3) 


3 X 10" ' 


15 


20,000 


0.9 


40 


0.02 (0.1) 


10-' 


50 



Table 1. The required number of stars for PSF cahbration of typical surveys, as predicted by our scaling relation 
(equation [5Q|) in function of the effective Signal-to- N oise Ratio of stars <Se f f (defin ed in eq nation [46|) and of the scientific 
requirement given in terms of (cTg^^)^ as defined bv lAmara fc Refregierl (j2007b) (i.e. the upper limit on the variance 
of systematics to ensure that systematics are sub-dominant compared to statistical errors when fitting cosmological 

S >arame ters) . The estimations of a[wo] and cr[wa] are made through the scaling relations proposed bv lAmara fc Refregier 
The mention 'nr' means that the constraint is 'not relevant'. Corresponding surveys (similar to the current 
generation of cosmic shear surveys) do not constraint the dark energy but rather constraint ag. However, (cTg^^ )^ 
does not depend on the considered cosmological parameter and thus these poor constraints on wq do contain all the 
information to compute (cTg^ )^. We consider that stars have an average Signal-to-Noise Ratio of 500 in space and 
300 from the ground. We also consider a fixed dilution factor (i^gai/^PSF)min of 1.5 for every survey. The bold line 
corresponds to the set-up of a DUNE-like survey (i.e. full-sky from space with medium depth and galaxy density) as 
taken as the reference in the scaling relation. Note that the surface of 20,000 deg^ (i.e. half of the sky) is the maximum 
surface of the sky relevant for cosmic shear. That is why it is often called 'full-sky'. In the same spirit, 15,000 deg^ is 
the maximum surface of the sky (relevant for cosmic shear) that can be covered from Mauna Kea with a reasonable 
air mass. 



parameter estimation. On scales smaller than the area con- 
taining A^* stars there is not enough information coming 
from the stars to calibrate the PSF. Therefore, the system- 
atic errors in the cosmological parameters may dominate 
over the statistical errors. This means that these small 
scales should not be used to constrain cosmology unless 
the variability is known to be extremely small. Our re- 
sults show that for current cosmic shear surveys this scale 
is about an arcminute, which may explain residual system- 
atics found on these scales in current analyses. In future 
all sky cosmic shear surveys, the data set will be increased 
by several orders of magnitude and a tight control of the 
PSF behaviour will be required to prevent this scale from 
being increased and to reach smaller scales. This places 
strong requirements on the hardware and observing strat- 
egy. For instance, for ground observations where the atmo- 
sphere prevents the PSF from being constant o yer several 
stars in a single exposure, Jarvis fc JainI ( 2004 ): Jain et al. 
(|2QQ6l ) suggest the shear correlation functions should be 
calculated by cross-correlating galaxies in different expo- 
sures. The PSF calibration error contribution would then 
be uncorrelated and averaged down. 

We also demonstrate how can depend on the com- 
plexity of the PSF. We define a 'complexity factor' for two 
different shapelet parametrisations and an elliptical gaus- 
sian parametrisation. This complexity factor is a number 
lower than the number of parameters used to describe the 
PSF and increases with the PSF complexity. There is a dif- 
ferent complexity factor associated with each PSF shape 
parameter. For accurate galaxy shear measurements, we 
find that the complexity factor for the PSF ellipticity is 
more important than that for the size. We summarise the 
dependence of iV* as a function of complexity, star Signal- 



to-Noise Ratio, galaxy size and cosmological requirements 
in a convenient scaling relation (see equation [5Q|) . Note 
that we work under the setup that the PSF can be per- 
fectly described by the model, thus we do not include in 
our error budget the systematics due to a poor choice of 
the PSF model. 

We find that pixelation degrades the PSF calibration 
accuracy. Our calculation of A^* holds for the optimistic 
case of infinitely small pixels (i.e. infinitely high resolu- 
tion) . A finite pixel size increases A^* . We have performed 
simulations to predict the increase in the case of an ellip- 
tical gaussian PSF. We show that for 2 pixels (or more) 
per FWHM, the rms scatter agrees with the standard de- 
viation predicted for infinitely small pixels to within 10%, 
which in turn can be translated into a 20% increase in the 
required number of calibration stars. For large pixels we 
see that the pixelation effect can be dramatic, for instance 
having 0.9 pixels per PSF FWHM would lead to a factor 
of 100 increase in the number of stars needed. This ef- 
fect may be mitigated to some extent by using dithering, 
which is beyond the scope of this paper. The main results 
of this paper are for the high resolution case that provides 
a strict lower limit on A^^^,. 
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Appendix A: Propagation of the PSF error 

In this appendix we detail how we derive equations [13] 
and [15] by propagating the error on each measurement of 
the PSF ellipticity and size to the error on the galaxy 
ellipticity and then to a^ys- 

The 2 component galaxy ellipticity €gai can be written 
in terms of radii R^^^ and Rp^^ and ellipticities eobs and 
€psF of the observed image and the PSF: 



jpohs ^PSF _|_ jpohs 



E^PSF 
^22 



€obs-^obs ~ ^PSF-^:^ 



PSF 



^obs 



RpSF 



(A.l) 



Consider now that we have estimators of epsF and Rp^p 
and the estimated values have small errors ^epsF and 
S {Rpsp) with respect to the true values. We can find the 
deviation from the truth of a single measurement of one 
ellipticity component (5egai,i by differentiation: 



de. 



gal,i 



dRp^p 



SRp^p 



depsF, 



-5e 



PSFA 



(A.2) 



which can be expanded to give equation [13] and when com- 
bined with equation [T2l 




^ (^psf) 



(A.3) 



which leads to equation [T5l when assuming that the ellip- 
ticity and size of the galaxy are uncorrelated (equation 
[T4|) and using equations [7] and [8] 
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-2.5 -2.0 -1.5 -1.0 -0.5 0.0 

Fig. 4. Non-gaussian PSFs used to study the effect of complexity in section [4?2l (lQ pixels per PSF FWHM, logarithmic 
color scale, the normalisation is such as the total flux is 10^). Each of them is described with a different shapelet model, 
described in table O 



label 




diamond 


Ncoef 


^^coef 


complexity factor 






(yes or no) 






^e(^K2) 


gauss 










a 


4 


no 


15 


2(3) 


3.27(2.58) 


b 


5 


no 


21 


2(3) 


3.27(2.58) 


c 


6 


no 


28 


3(4) 


4.47(3.74) 


d 


7 


no 


36 


3(4) 


4.47(3.74) 


e 


8 


no 


45 


4(5) 


5.66(4.90) 


f 


9 


no 


55 


4(5) 


5.66(4.90) 


g 


10 


no 


66 


5(6) 


6.83(6.06) 


h 


11 


no 


78 


5(6) 


6.83(6.06) 


i 


12 


no 


91 


6(7) 


8.00(7.21) 


j 


4 


yes 


9 


1(3) 


1.63(2.58) 


k 


6 


yes 


16 


2(4) 


2.83(3.74) 


1 


8 


yes 


25 


3(5) 


4.00(4.90) 


m 


10 


yes 


36 


4(6) 


5.16(6.06) 


n 


12 


yes 


49 


5(7) 


6.32(7.21) 


o 


14 


yes 


64 


6(8) 


7.48(8.37) 


P 


16 


yes 


81 


7(9) 


8.64(9.52) 



Table 2. Properties of the PSF models investigated in section [4^ The first one is the elliptical gaussian (labelled 
'gauss'). For each of the 16 other basis sets, we study one PSF (lab elled a to p) shown i n figur e [H From left to right, 
the columns are: nmax and diamond defined in section [331 and bv iMassev fc Refregierl (|2005l ): the number Ncoef of 
coefficients in the basis; the numbers A/"^^^^^^ of coefficients with m = 2 (and m = 0) in the basis and the complexity 
factors defined by equations [39] to [411 



Appendix B: Complexity factors for shapelet of complex coefficients /n,m (see equation [38]). The scale 

models parameter P and the center of the basis can be tuned to 

get a sparse description of the object. Most of the time. 
As stated in section [331 any object can be described by a good choice is to take /3'^ = and the center of the 

the sum of polar shapelet functions Xn,m weighted by a set basis corresponding to the centroid of the surface bright- 
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ness. Then, for simple objects like stars with reasonable 1. We neglect the term proportinal to in equations IB. 91 
substructures and tails, the required number of coefficients and lB.lOl This is justified by the 'small ellipticity' as- 
is typically less than 20. The 'diamond' configuration, de- sumption adopted all along this paper, 
fined for nmax even, consists in setting fn,m = when 2. We choose = 2/3'^ in equations \B.6\ [B?7l and [R8l 
\m\ > nmax - n. 

For given values of and the center of the bas is, equa- 
tions 50, 54 and 55 of iMassev fc Refregier (|2005h can be 
summarised by: 



= VI^p X So 



= 2/32 X Si/ So 

= 82/81 

with: 

''T'max 

50 = ^ fnO 

n even=0 

51 = £ (n + l)/„o 

n even=0 

''T'max 

52 = \/^(« + 2)/„2 . 

n even=2 

From equations IB. H [R2l and [2H it follows that: 
1 



a[R^] = -x2P^x^j,2 



a[e] 
with: 

= 



1 ^ 



N/2 



n even=0 

and, for the non-diamond configuration: 



(B.l) 
(B.2) 

(B.3) 
(B.4) 
(B.5) 

(B.6) 
(B.7) 

(B.8) 



^ £ [n{n + 2) + eHn+lf] 



n even=0 



N/2 



N{N + 4:) 2fN,^, , 
^ ^ ( Y(^ + 4) + 1 



(B.9) 



or alternatively, if the diamond configuration is used: 



nmax -2 nmax 

^ ^ [n(n + 2)]+e^ ^ [(n + 1)^] 



7V/2 



^-(^--2) 



n even=0 



n even=0 



■ I y (Af + 4) + 1 



(B.IO) 



where N is the largest even integer lower than or equal to 
^max- Note that a[e] is the standard deviation of one com- 
ponent of €, according to our definition given in equation 

1 

To obtain equations [39] to [411 we make two simplifica- 
tions: 



